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Abstract. - The influence of an initial momentum on the appearance of "quantum resonances" in 
the delta-kicked rotor system is explored experimentally. We show that for certain initial momenta, 
a resonance can be negated entirely, whereas at others a resonance can be made to appear. At a 
larger number of kicks, all resonances are shown to narrow. We show that as a consequence, the 
individual "diffraction peaks" may split. We compare our results to numerical simulations as well 
as analytical theory. 



The delta-kicked rotor system has been a paradigm for the study of "quantum chaos" : 
quantum systems for which the classical equivalent exhibits chaos. A model system that 
has been studied extensively is the "atom optics" kicked rotor, which consists of a two-level 
atom in a detuned, standing wave of laser light. This system has been studied extensively 
theoretically [1] and experimentally [2,3] for cold atoms originating from a Magneto-Optical 
Trap (MOT) . The temperature of the atoms in a MOT is in the order of microKelvins, much 
larger than the heat input from the recoil of a single photon, and the typical dimensions of 
a MOT is in the order of a millimetre, much larger than the wavelength of the light in the 
standing wave. Experiments have demonstrated dynamical localisation [4] and the existence 
of "quantum resonances" [5], depending on the kick period. 

At these quantum resonances, an effect similar to the optical Talbot effect [6, 7] takes 
place. In this effect, the near field after a light field diffracts from a grating is recreated at 
a certain distance away from the grating. The analogue in atom optics is observed when 
a diffracted plane deBroglie wave is diffracted after a certain "Talbot" time Tt = it/(2lur) 

the near field is recreated. Here lor — is the recoil frequency, and we have defined the 
laser wave number = 2tt/\ and the mass of the atom m. 

When the kick period T corresponds to an integer times the Talbot time (T — tjiTt, 
with integer m), the accumulated phase of the wave function due to moving with a velocity 
2h~kL is exactly 2ir, and the kicks add constructively. This leads to a quadratic growth of the 
variance of the velocity distribution, or energy, with the number of kicks. For half integer 
times the Talbot time (T = (m + |)2t), The kicks add destructively, and the effect of the 
first kick is canceled by the second. However, for velocity distributions much broader than 
a photon recoil, an average over the initial velocities leads to an observed resonance with a 
linear energy growth for both integer and half integer times the Talbot time. 

Recently, experiments [8] and theory have turned to using atoms originating from a Bose- 
Einstein Condensate (BEC), which are orders of magnitude colder, less than the heat input of 
a single photon. Experiments have shown that for these conditions, higher order resonances 
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can be observed [9], which have been shown to depend sensitively on the initial momentum 
of the atoms. Furthermore, a range of experiments have exploited ratchets [10,11], and most 
recently, the dependence of its operation on the initial momentum [12, 13]. 

Here, we present the first measurements of the effect of the initial momentum on the 
primary quantum resonances in the delta-kicked rotor. We use an ensemble of atoms with 
a momentum spread much less than a single photon recoil, which is kicked at kick periods 
corresponding to the quantum resonances. We explore the first three quantum resonances, 
for two and four kicks. Furthermore, we investigate details of the velocity distribution after 
a series of kicks. 

To describe our system, we are considering a quantum particle of mass m moving in 
one dimension in a pulsed, spatially periodic potential with period A/2. The Hamiltonian 
is given by 

2 N 

H = — + h<j) d cos 2k L xV^ 5{t-nT) (1) 
2m ' 

n— 1 

where T is the pulse period and <pd is the kick strength. It is well known that for a quantum 
system, strong resonances exist for the Talbot time T = ttlTt where the effects of the kicks 
add constructively, and strong diffraction occurs. Furthermore, at T = (m + i)7r/(2aj^) 
anti-resonances occur, where the effects of the subsequent kicks add destructively, and only 
very little diffraction occurs. In the experiment, the delta-function is represented by a block 
function with a width r <C T and an area (pa- 

We can describe the time evolution of the quantum particle during one kick and the 
subsequent free evolution period with a Floquet operator: 

p _ e -ip 2 T/2mft e -i<Pd cos(2k L x)r /h (2) 
v v 

= Ufree = ^kick 

We simulate the experiments using the well-known split operator method for the Floquet 
operator, which is ideally suited to this experiment, with 2 16 grid points. We start the 
simulation with a minimum uncertainty wave packet, with an initial momentum Tiki. The 
initial wave function is given by 

1 - - 2 

V-(x) = e ^l e ~ ikiX (3) 



2tt<j u 

The width a w — 5 /an of the wave packet is chosen so that it covers several wavelengths of 
the kick laser. The simulation is run with a distribution of initial momenta, corresponding 
to the initial momentum distribution found in the experiment. These results are added 
incoherently to obtain a final momentum distribution. 

We can condense the effects of the kicks on the velocity to a single number, the variance. 
This is proportional to the average kinetic energy of the atoms, and can be expressed in 
recoil energies. The energy is not influenced much by the details of the original momentum 
distribution, as the typical energies are larger than one recoil, whereas the initial energy is 
much smaller than one recoil. 

For kick periods that are a half integer times the Talbot time, an analytical expression 
has been derived for the amplitudes of the momentum states with momentum 2jhk after n 
kicks by [14] and [15]: 

Cj = Jj (^^taT^) i j e- ij ^ r e-^ 1 (4) 

with T = + 2/3)1, (3 the quasimomentum in units of 2KkL, and integer I = 2T/Tt- 
From equation 4, the second order moment of the momentum distribution can be found: 

(p 2 ) = (2ft**) 2 £ \ Cj \ 2 (j + P) 2 (5) 
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Fig. 1: Absorption images after two kicks with a kick period of T = Tt/2 for different initial 
momenta in units of photon recoils. 

The values obtained from equation 4 are in excellent agreement with those obtained from 
the simulation, and with the experimental data. 

In our experiment, we use ultracold atoms sourced from an all-optical BEC of 87 Rb 
atoms at a temperature of 50 nK, which corresponds to a velocity distribution with a a w 
0.25 recoils. The BEC is formed in a crossed pair of CO2 laser beams. The ~ 5 • 10 4 atoms 
in the BEC are in the F = 1 hyperfine state, and spread over all magnetic substates. The 
apparatus is described in detail in [16]. A pair of counter-propagating diode laser beams, 
overlapped with the BEC, creates the standing wave potential. In order to reduce mean 
field effects, the atoms are left free to expand by extinguishing the CO2 laser for 500 /is 
after forming the BEC before starting the kick sequence. 

The horizontal grating laser beams are derived from a home-built diode laser system. 
The laser is locked to the F = 2 -> F' multiplet in 85 Rb at 780 nm, and thus detuned by 
2.2 GHz from the relevant F = 1 — > F = 2 transition in 87 Rb. The laser beam is split 
50%, and both parts are passed through a separate acousto-optic modulator (AOM). The 
AOMs are driven by an Arbitrary Function Generator (Tektronix AFG-3252), amplified by 
home-built RF amplifiers. The AOMs are pulsed on simultaneously, but with a tunable 
frequency difference Aw. This frequency difference gives us an effective initial momentum 
Pi for the atoms, given by Pi/p re c = Aw/(4wr). It should be noted that what is important 
is the relative phase of the sine waves driving the AOMs in the subsequent pulses, caused 
by the small frequency difference, not the actual frequency difference. 

Care is taken that the laser pulses arrive at the experiment at the same time. A shutter 
is used in the laser beam to ensure that it is totally extinguished during the evaporation 
phase to create the BEC. The lasers beams are passed through polarisation preserving single 
mode fibres to the experiment. The total laser power at the experiment can be varied up to 
2 mW in each beam. 

The laser beams are focused to a 100 /im spot at the BEC, in a standing wave config- 
uration. The foci are much larger than the size of the BEC (<~ 10/xm), thereby yielding a 
constant interaction strength over the BEC, whilst increasing the intensity. The increased 
intensity allows us to perform the experiment at large detunings to suppress spontaneous 
emission. 

An absorption image of the atoms after a 5 ms expansion time is taken to obtain the 
momentum distributions. The 2D momentum distributions are summed over the width of 
the cloud to obtain a ID momentum distribution. In all cases, the experiment is repeated 
three times. 

In figure 1 we show absorption images of the atoms after two kicks at the antiresonancc 
T = Tt/2 and after an expansion time of 5 ms for different initial momenta, as set by the 
frequency difference in the kick beams. At p — we see almost no diffraction, as expected 
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Fig. 2: The number of atoms as a function of final momentum, for a range of initial momenta from 
to 2hk. Two kicks with period of Tt/2 (left) and T = Tt (right) are applied, with experimental 
profiles (dotted,red)and results from the simulation (full line) for 4>d ~ 1 

from the Talbot effect. At pi = 1 p YCC , as well as at pi = 3 p rec , we see strong diffraction, as 
Tt/2 corresponds to a resonance for these atoms. At pi = 2 p rcc , as well as at pi = 4 p lec , 
there is again no diffraction. 

Next, we vary the initial momentum pi in 8 steps from to 2 p rcc . On the left in figure 2 
we show the resulting ID momentum distributions, averaged over the three repeats of the 
experiment. It can be seen that at zero velocity, there is indeed an anti-resonance and very 
little diffraction occurs. However, it turns into a resonance at pi = p rcc , with significant 
diffraction, and back to an anti-resonance at pi = 2p lec with little diffraction. Also shown 
in the figure are the results of the simulation. When experiments are performed with a 
broad initial momentum distribution, an average over many different Uj yields an observed 
resonance at T = Tt/2, as has been shown in many publications, see, e.g., reference [5]. 

On the right in figure 2 we show the resulting ID momentum distributions for two kicks 
at a spacing of T — Tt, for varying pi, both experiment and simulation. At pi = 0, there is a 
resonance as expected, with significant diffraction. However, at pi = 0.5p rec , this turns into 
an anti-resonance, back to a resonance at pi = p rec , back to an antiresonance at pi = 1.5p rcc 
and finally to a resonance at pi — 2p rcc . Again, averaging over a range of initial velocities 
will show the same resonance at T — Tt, even though the dependence on the initial velocity 
cycles at twice the rate. At T = 3Tt/2, the cycle of the amount of diffraction varying with 
the initial velocity is at three times the rate we see for T = Tt/2, and so on. The simulated 
curves show good agreement in terms of the relative heights of all the diffraction peaks. 

A characteristic feature of these velocity distributions is the variance of the momentum 
distribution, or the mean kinetic energy in units of the recoil energy E r — Tkjfj. In figure 3 
we show the variation of the energy on the initial momentum for T = Tt/2 (a), T = Tt (b) 
and T = 3Tt/2 (c), experiment and simulation. For T — 3Tr/2, the energy varies at three 
times the rate for T — Tt/2. For the data points, the averaged profile is fitted to a number 
of Gaussians, one for each diffraction order, to obtain the energies. The uncertainties are 
generated by fitting each experiment to obtain the energy and taking the standard deviation. 
It was found that a direct numerical variance of the velocity distributions gave similar results, 
but much less consistent due to its sensitivity to small noise peaks at large momenta, and 
consequently these values were not used. 

For a larger number of kicks, it is well known that the width of the quantum resonance, 
when plotted versus the kick period, gets smaller as the number of kicks gets larger. Expect- 
ing a similar behavior, we investigate the dependence of the energy on the initial velocity 
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Fig. 3: On the left, the final energy as a function of the initial momentum for two kicks at a pulse 
period of 33.15 /xs (a), 66.3 /xs (b) and 99.45 /is (c). On the right, the energy after four kicks at a 
period of Tt/2 (d) and Tt (e). The experimental data and uncertainties (markers) are shown as 
well as the results from the simulation (line). The kick strength ^ fs 1 for these figures, except 
that <f> d = 0.9 for (c) 

for four kicks. In figure 3 we show the energy as a function of the initial velocity for four 
kicks, for both the anti-resonance (d) and the resonance (e). For T — Tt/2 (d) there is an 
anti-resonance at zero momentum, yielding small energies. There is a small maximum at an 
initial momentum of j recoils, decreasing again to close to zero at ^ recoils. At an initial 
momentum of one recoil, there is a stronger maximum in the energy, decaying again to small 
energies at | before another smaller maximum and returning to small energies at two recoils 
initial momentum. All these features, prominent in our calculations, are reproduced by the 
experiment. 

For T — Tt (e) a larger number of oscillations is observed in the simulated curve, along 
with more narrow maxima. The experiment reproduces the maxima, but at current cannot 
resolve the small-period oscillations for these settings. 

Now, we will take a closer look at the momentum distribution after four kicks, with 
an initial momentum of one recoil. The initial momentum distribution from the BEC has 
a certain width a = 0.18 recoils, which qualitatively alters the momentum distribution 
after the kick sequence. Ideally, the initial momentum distribution is a delta function at 
Pi. In reality, the initial momentum distribution has a width, and there are also some 
atoms at a momentum pi + Ap, where Ap is much smaller than one recoil. To illustrate 
this, we performed a simulation for a range of initial momenta. In figure 4 we show the 
probability (grayscale) of finding the atom at a final momentum (horizontal axis) for a 
range of initial momenta for both the antiresonance and the resonance. We also display the 
initial momentum distribution of the BEC (black line). The sensitivity of the system to the 
initial momentum translates into peaks coming into existence at a small offset momentum 
Pi + Ap that are not there at pi (pi is equal to one recoil in this case), or, conversely, a peak 
disappearing at pi + Ap which was present at pi. 

As mentioned earlier, we model the effect of a finite width of the initial momentum 
distribution by carrying out the simulation for a range of Ap 7 and we subsequently add 
the momentum distributions thus obtained incoherently, weighted by the height of the mo- 
mentum distribution of the original BEC at Ap. In effect, to obtain the final momentum 
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Fig. 4: The probability of finding atoms (grayscale) at a final momentum (horizontal axis) for a 
range of initial momenta (vertical axis). Also shown is a Gaussian momentum distribution which 
represents the distribution of the initial BEC. This is for four kicks at T = Tt/2 (top) and T = 
Tt (bottom), average initial momentum (pi) = p roc and d = 1 

distribution, we take horizontal profiles in figure 4, and sum these, weighed by the initial 
momentum distribution. 

These momentum distributions are shown in figure 5. The difference in width between 
the different "diffraction orders" can easily be observed and related to the separation of one 
velocity peak into multiple peaks in the simulation. For instance, the peak at a momentum 
p = 1 recoil is split into three, with varying heights in the two situations. This would be 
caused by the peak disappearing for Ap — 0.2 recoils, and subsequently re-appearing for 
Ap = 0.4. The latter peaks appear in the final momentum distribution even though they 
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Fig. 5: The momentum distributions after four kicks for T = Tt/2 (top) and T = Tt (bottom) with 
the experimental profiles (full lines) and the simulation (dashed lines). In both cases, (pi) = p ICC - 
Again, 4> d « 1 
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are not strongly weighted by the initial velocity distribution. 

The peaks at p = (3, —1) recoils are split into two, which would be caused by peaks 
appearing at a small Ap. On the other hand, the peaks at p = —(5, 7) recoils are considerably 
narrower than the initial momentum distribution, which is caused by the fact that this 
diffraction peak only appears for very small Ap. These observations can be verified in 
figure 4, by drawing a horizontal line at fixed Ap and finding the diffraction orders. 

At a larger number of kicks, the simulations show that the width of resonances as a 
function of initial momentum will get smaller. Determining the energy after a number of 
kicks could be used to measure the velocity distributions of an atomic sample with great 
accuracy. As it is influenced by the velocity only, it is not influenced by the initial cloud 
size as is the time-of-flight method. 

In summary, we have for the first time experimentally demonstrated the dependence 
of the quantum resonances in the delta-kicked rotor on the initial velocity of the atoms. 
We show a sinusoidal dependence of the energy on the initial momentum for two kicks, 
and a more complex behavior with the same period at four kicks. We propose to use the 
momentum dependence of the kicked rotor to select narrow parts of an initial momentum 
distribution, or to measure a given momentum distribution with great accuracy. 

The authors would like to acknowledge fruitful discussions with Scott Parkins regarding 
this work. 
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